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Purpose: To relate the physical transmittance parameters of the water equivalent 
Gafchromic EBT 2 film with the delivered dose in a transparent absolute calibration 
protocol. The protocol should be easy to understand, easy to perform, and should 
be able to predict the residual dose error. 

Methods: The protocol uses at most two uncut calibration films for a whole 
batch. The calibration films are irradiated on different positions with static fields 
with varying dose levels. 

The transmittance trough the film, T, is preferred over the optical density. An 
analytical function is used to associate T with the delivered dose based on physical 
characteristics of the film: the minimal and maximal T- values, To and T^, and a 
factor scaling the dose, ^3. 

The dose uncertainty of the protocol is calculated using an error propagation anal- 
ysis of the calibration curve. This analysis makes use of the small scale and the large 
scale variations in the film response (local and global uniformity). Before calibrat- 
ing the film, the stability of unirradiated films is evaluated. The relevance of the 
positions of the different dose levels is studied. 

Results: Both a production and a spatial dependency of the unirradiated films 
was noticed. This spatial dependency required an avoidance of a quarter of the film. 
The transmittance is also affected by the storage and the transport conditions. 

The lowest residual errors and the highest significance for T , T^, and (3% was 
achieved with arbitrary dose levels for different positions on the film. Two calibration 
films give the better results. Significant calibration curves were found for both the 
red and the green color channel. 

Large differences (> 0.02Gy) between calibration curves, with different positions 
for the dose levels are seen. 

In the [0.04,2.5]Gy dose range the red channel calibration curve has dose errors 
ranging from -2.3% to 4.9%. For the green channel the error range is [-3. 6%, 6. 9%]. 

The calculated error propagation was able to predict an upper limit for the red 
channel dose errors. For the green channel this prediction was not achieved. 

Conclusions: The gafchromic EBT2 films are properly calibrated with an ac- 
cessible robust calibration protocol. The protocol largely deals with the uniformity 
problems of the film. The proposed method allowed to relate the dose with the red 



channel transmittance using only To, T^, and a dose scaling factor. Based on the lo- 
cal and global uniformity the red channel dose errors could be predicted to be smaller 
than 5%. 
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I. INTRODUCTION 

Radiochromic films, such as Gafchromic EBT and EBT 2 films^, are a popular tool to 
measure dose depositions from ionizing radiation. Most of the previously published studies^ 
have focused on the gafchromic EBT film, the predecessor of the EBT 2 film. However, the 
active layer of both types of films is the same, which results in comparable outcomes for the 
EBT and the EBT 2 filmP^l. The main advantages of these type of films are the fact that 
these films have the following properties: 

a. A near flat response to different photon energies^ 

b. Close to water equivalence (Z eff EB t 2=6. 98 4 ) 

c. A thickness less then 0.3mmP31 

d. A low perturbation of the irradiated media, which is certainly important for the 
evaluation of recent rotating treatment techniques like VMAT 8 

e. A reasonably large dynamic range ([0,10]Gy 6 , [0.1,8]Gy^) 
/. A high resolution (0.22/zmPEl) 

Disadvantages of the films are their non-uniformity and production stability, which has 
an impact on the choice between transmittance (T), and optical density (OD), as dose 
dependent quantity (section [II A 1 ). 



However, the combination of the afore mentioned properties makes gafchromic EBT 2 
film, further referred to as film, usable in almost every application in radiotherapy. For ex- 
ample, at the University Hospital Leuven the films are clinically used : for in-vivo measure- 
ments of total body and electron skin irradiations; for Volumetric Modulated Arc Therapy 
quality assurance; for internal audits of brachytherapy, and stereotactic treatments; for the 
introduction of new treatment techniques^, and new clinical trials (FLAME triaP^). 

Additionally, there is a large range of research purposes such as; out of field dosimetry^, 
estimation of microscopic dose distributions in nano-particle enhanced radiation therapy^^, 
and the optimization of the dosimetry of fluoroscopic CT protocols. 

The wide use of these film requires a transparent and accessible absolute calibration proto- 
col. The calibration methodology aims to use one or two uncut calibration films to represent 
a whole batch. The use of an uncut film avoids puzzling with calibration film fragments, 
and avoids sharp transmission and diffraction edges in the film scanner transmission system. 
Calibration films and patient quality assurance (QA) films are always scanned successively, 



but they are not necessarily irradiated on the same day. An error propagation analysis of 
the calibration protocol is built using the small scale and the large scale variations in the 
film response. 

II. MATERIALS AND METHODS 

The calibration protocol is based on eight static fields irradiated on an uncut 8x10 inch 
film, see figure [Tj Table [I] lists all the films used in this work. 

To relate T with the dose, D, only three parameters are used (/3i, (3 2 , fiz). These parame- 
ters can be reduced to the zero dose transmittance, T , the infinite dose transmittance, T^, 
and a factor scaling the dose, 4-. 

A. Theory 

1. Optical density versus transmittance 

Historically, OD is the quantity of choice because of its linear relation with the dose 
delivered, e.g. TG 55 reports a linear net OD response from to 30 Gy, and the from 30 to 
lOOGy for the MD-55-2 radiochromic filmP^. Such a linear relation allows a scaling of relative 
dose measurements. For the gafchromic EBT film, Rink, Vitkin, and Jaffray] mentioned a 
nonlinear OD response^. 

An exponential function can be used to describe such a nonlinear relation between the 
dose and OD^"^ , e.g. OD= A ■ (l — e~ BD ). Such an expression results in zero optical 
density, when no dose is applied, OD 0Gy = 0. On the other hand, the optical density is 
defined as OD= log 10 ( t 1 ) (f^, formula 12.34). Where <pi is the visible light fluence incident 
on the film and (f>t is the fluence transmitted through the film. To reconcile both expressions 
it is necessary to replace fa by the fluence transmitted through an unirradiated film, 0*, 
which results in zero OD when no dose is delivered (OD G y = 0). 

Over the course of 2010 we had variable success using the scan of a separate unirradiated 
film as an estimate of <fi* (table ffl Film^-B^). Using this approach two main problems, 
were noticed. Firstly, a stability problem; the unirradiated film, and therefore the estimated 
</>*-value, changes over (production)-time. And secondly <fi* is spatial-dependent, see figure 
[2} This non-uniformity requires a location dependent 0*-estimation. Because a stable 



estimation method for 0* was not obvious, and because the loss of linearity for the OD-Dose 
response, our research concentrates on a straightforward association between the physical 
quantities transmittance (T = &) and dose. 

2. Association between transmittance and delivered dose 



Figure [3] illustrates the association between the dose delivered and the transmittance 
as measured by us. Mathematically, T decreases monotonically as D increases. Such an 
association can be described by a rational functional, e.g. T = g+ D ■ This function is 
more interesting than a polynomial function due to the absence of points of inflection. The 
rational function is extended with a term depending on the pixel location, the /Jj ■ X-term 



in equation la This term takes into account a non-uniformity along the 10 inch side of the 
film. The transitions from equation la to lb can be found in appendix |A~} 
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Ignoring the non-uniformity term, (3^ = 0, equation [la] can be interpreted as follows 
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The transmittance when no dose is delivered equals the ratio of (3i and f3%. 
Too = 02 '■ 02 is the transmittance resulting from a theoretical infinite dose, T^ : 



^ , M^ 
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(l'Hopital's rule). In appendix 
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equations 



la 



and 



lb 



are rewritten to 



equations 



2a and |2b[ illustrating the relation of the physical quantities T, T , T^, and D. 
From these equations, J- can be seen as a factor scaling the impact of the dose. 
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3. Error analysis 



The partial derivatives of equation lb are used in an error propagation analysis. This 
analysis estimates the error made when film measurements are performed, using the proposed 



calibration methodology. The dose uncertainty of the protocol, an, is calculated using 
estimations of the position variance ax, the transmittance variance Ct, and their covariance 



Cov (T,X) (equation 3a 3d). 
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B. Film handling 



1 . Digitization 



The films were positioned in the middle of an Epson 10000XL flatbed scanner, the 8 inch 
side is parallel to the long edge of the scanner bed. Uniformity artifacts due to the position 
on the scanner are well knowrp2ZE therefore, a plastic positioning tool ensures the same 
position for each scanned film. The films are studied with the orientation mark, clipped by 
the manufacturer, directed to the lower left corner of the scanner (figure [I]). 

To measure the visible light fluence, <p, the scanner is equipped with a transparency cover. 
Scanning is performed using the Epson scan software, with all corrections turned off. A fixed 
8x10 inch area is scanned in 48-bit color mode with a color depth of 16-bit; ie. the pixel 
values range from to 2 16 — 1. Files are saved in tagged image file format (TIFF) with 150 
dpi resolution. These TIFF images have three layers containing red, green, and blue scan 
values. All subsequent film processing is performed using home made Matlab scripts (The 
MathWorks, Inc). 

As described in the previous section, transmittance values are determined by the ratio of 
the fluence transmitted through the film and the fluence incident on the film Trgb ~ " ; '" :Ii 



'iRGB 



The definition is valid for all three color layers, this is indicated by the RGB index. 0iRGB is 



measured by a scan without a film the resulting signal is within 1% of the maximal fluence 
value (2 16 — 1). Therefore, 0jrgb is set constant to this value, and : 

Trgb = ^ RGB (4) 

Our films are scanned at least 48 hours post irradiation, which is the preferred recom- 
mendation in the TG 55 reporfPJ Andres et alP reported variations in OD smaller than 3% 
between 2 and 24 hours postirradiation. They report this value for irradiations up to 2.5Gy. 
Taking in to account that higher doses have a slower density growthP, and the decrease of 
the subsequent density growth as a function of time, 48 hours between irradiation and read 
out guaranties differences in density growth below 3%. 

2. Film storage 

The film response is affected by the storage environment, e.g. the impact of moisture 
and temperature on optical density 22 . Therefore, all films are separately stored in folders 
of a file cabinet, which is placed in an air conditioned environment. The file cabinet limits 
exposure to daylight. The separated folders and air conditioning ensure identical exposure 
to environmental conditions. The relative moisture and the temperature of the storage room 
are monitored over the course of a week. 

A set of fourteen films was used to assess the effect of the moisture and temperature 
conditions of the storage environment. The manufacturer supplied the films, stored in 
vacuum packing. The vacuum pack was chosen to exclude environmental effects during 
shipping and storage by our local distributer (PEO, Radiation Technology, Hoogstraten, 
Belgium). The first seven films (Film^-vv) were separated by an interleaving tissue, the 
last seven were not (Filmy 7 _y 14 ). The vacuum storage was broken on October 20 th , 2010. 
The films were scanned unirradiated on the same day and on subsequent days to evaluate 
the environmental impact. 

For the evaluation of the storage effect, local transmittance variations (variations on 
short distances - noise) were not addressed. Therefore, the transmittance layers is sub 
sampled by taking the average over 5x5 pixel blocks^. This operation results in the RGB- 
transmittance layers denoted by T^ GB . The effect of our local storage conditions is then 
evaluated by calculating the average of T BGB for the different films, Film^-v^- The outer 
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20 pixels were excluded from the calculation to avoid effects of the marks as well as other 
effects at the film edges. We denote the average and corresponding standard deviation as: 

Orgb)* ± SD (T£ GB )* (5) 

The (.}* and SD (.)* notation is used to indicate the difference between the average of a 
whole film, and the average of a small regions of interest (ROI), denoted by (.) and SD (.) 
(section II C 2 , equation [6]) . 



C. A static field calibration methodology 

1. The calibration field 

The calibration films are irradiated with 4x4 cm fields with a 6MV photon beam from 
a Varian Linac 2100C/D. Figure [l] shows a schematic overview of a calibration film. The 
sizes, the location and orientation of the calibration fields are chosen to maximize the sepa- 
ration between fields. The selected field size (FS) guaranties lateral electronic equilibrium. 
According to Todorovic et a/.^the calibration curve is not affected by the field size, for field 



sizes ranging from 2x2 to 10x10 cm. 

The protocol uses static fields, to simplify the physics of the calibration fields. For 
example, by avoiding the use of a multileaf collimated (MLC) field, leaf transmission effects 
are excluded. However, the irradiation of a single calibration field still generates a low dose 
contribution to the positions of the other fields due to scatter and/or leakage, see section 
IHC31 

2. Transmittance ((T RGB )) 

The calibration film is positioned in a plastic water phantom. The phantom is a stack of 
twenty 40 x 40 x lcm RW3 plates, which is positioned with 95 cm source to surface distance. 
The film was positioned at 100cm from the source. This way there is sufficient build up and 
backscatter (5-15 cm). 

On the calibration film, eight calibration segments are defined. These eight segments 
are 1 x lcm ROI's situated in the center of the 4 x 4cm diamonds as shown in figure [1} 
The segments' edges are parallel to those of the diamonds. The segments can be irradiated 



with a known dose or can be left unirradiated. The average T-values of the segments are 
calculated; 

(T RGB ) ± SD (T RGB ) (6) 

The different calibration segments are irradiated by applying couch shifts. The monitor 
units (MU) are chosen to have an equal spread between and 100MU, and between 100 and 
225MU, see table |T] The dose range of the calibration curve is ~ 0.4 to 2.5 Gy, which covers 
our standard external beam dose range. 

3. Leakage and scatter conditions (Di ic) 

In order to account for leakage and scatter contributions to the individual calibration 
segments, we use the following methodology and assumptions: 

a. The different segments have symmetric contributions i.e. scatter and leakage contri- 
bution of segments 1 to segment 8 is identical to the contribution of segment 8 to segment 
1, apart from a scaling factor depending on the delivered dose. 

b. The ionization chamber (IC) response is weakly dependent on the energy spectrum, 
therefore, IC-measurements are considered to be accurate enough for this purpose. 

For IC measurements, a farmer type IC (FC 65-G TNC, SN 752) is combined with a 
SI (Standard Imaging, Inc., Middleton, WI USA) electrometer (SI SN 070112). The IC is 
placed in the RW3-phantom, located at the center of the measuring plane, at 100cm from 
the source. 

The eight dose contributions to segment i, D it j ioomu, are measured by placing the IC 
in the center of segment i while irradiating 100 MU to each segment j, j 6 [1,8]. Again, 
couch shifts are used to position the IC in the center of the different segments. The eight 
measurements should only be performed once, after which only the dose of a reference field 
should be measured, D iiCal 100 mu- 

This reference measurement is repeated before the irradiation of a calibration film. Sub- 
sequently, the dose delivered to segment i, A ic (equation [71), is the sum of the contribu- 
tions of the eight segments on segment i, scaled by the ratio of the output measurements 

/-Dj.iCal 100 MU \ 



Di : i 100 MU ' 

c 
n A.iCal 100 MU >T^ MUj 

Di IC " a^oomu ' \ mm x D « 10 ° MU (7) 

3 =l 
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Where MUj denotes the MU delivered to segment j. 

D. The calibration curve 
1. Location dependence 

Recently Micke, Lewis, and Yu reported a lateral-artefactP', which is a location depen- 



dence of the calibration curve of gafchromic films, ie. the /^-parameters of the calibration 
curve are affected by variations of the different dose level positions. To study the impact of 
such an effect on our calibration methodology, eleven calibration films are irradiated with 
varying dose levels for the different segment locations. This is done by altering the order of 
MU used to irradiate the calibration films (Filing __£ n , table [I]). 

Increasing MUs with increasing segment number are used thrice on two different days, 
Film^ and FimiE 2 £ 3 . Decreasing MUs with increasing segment number is used for Filing. 
Filing has an arbitrary spread of the MUs on the different locations. Filing has the same 
MUs as Filing but an arbitrarily segment is left blank and segment six got 225MUs. Film^ 
has a low and a high dose field on each row. 

2. Calibration parameters estimation (Phigb) 

Subsequently, the (T RGB )— values, and the corresponding Z) IC — values, of these eleven 
films are pooled in ten data sets. 

Three pooling methodologies are considered. The first method (data set 1), pools the 
segments of all eleven films. The second method (data sets 2-8), pools the segments of a 
single film: respectively Filing, Filing, . . ., Film^, and Film^ 10 . The third methodology 
(data sets 9,10), pools the segments of two films, Filing + Filing, and Filing + Filing. 

The two film combinations are introduced to increase the number of data points. The 
Filni£; 2 + Filing-combination only enlarges the number of data points, while the Filing + 
Film^-combination ensures both an enlargement of the data set as a varying distribution 
of low and high dose values on the different segment locations. 



For each data set, the /3jRGB-parameters (equation la) are estimated for all three color 



channels. The estimation is performed in R (www.r-project.org) using a non linear least 
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square fit of equation [Tal The t-statistic and the corresponding two sided p-values are 



calculated to evaluate the significance of the /3jRGB-estimations (a = 0.05). 

3. Six and Eight segment approach 

Because of film non-uniformities (Figure [2]) all preceding methods were repeated including 
only six segments. The segments with X = —9cm are avoided in this six segment approach. 
The methods described before, including all segments, will be referred to as eight segment 
approach. Only included segments are used in the further evaluation of the approaches (e.g. 
residual dose error calculation etc.). 

E. Error analysis 

1. Residual dose error, e 

The /3jRGB-estimations of the previous section are used to convert the (Trgb) — values to 
dose (equation [81). To compare the eight and the six segment approach the residual errors 
are calculated (equation [9]). 

(& • (T) -03-P*-X)-j3 1 



D Film rgb p 2 _ {T )+p 4 .X (8) 

x ' ' RGB 

^RGB = -DFilm RGB — -Die (9) 

The six segment approach on the Filing + Filing data set has both the lower residuals 
as the lower p-values for the /3jRGB-estimations. Subsequently, the calibration curve resulting 
from the Film^ + Filing data set is used to calculate the dose values for all the segments of 
all the films, Filing _E n . The resulting residual dose errors, £rg, are qualitatively compared 



with an estimation of the dosimetric variation, op 3a The blue color channel is not evaluated 



because of the non-significance for the /SjB-estimations. 



2. Dosimetric accuracy (&d) 



As described in equation [3a] the dosimetric accuracy consist of components involving 
positional accuracy (crx), local uniformity or transmittance accuracy (<tt), and finally a term 
taking into account variations of T as function of the X-position (Cov (T, X)). No significant 
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variations of the transmittance as function of the F-position were noticed (ANOVA, a = 

0.05). The contributions of the different terms are determined as follows; 

ax: A 1mm positioning error is assumed both during irradiation as scanning. These 

two positioning errors are considered to be independent. 

Or-' The local uniformity is defined as the relative standard deviation of TrgePI 
( r.gb) _ aT rgb j g es ti ma ted by the maximum of the standard deviations of all the 

segments of all films, Film^^^: 

(J T RG B= p .,niax (SD(T RGB )) (10) 

Cov (T, X) : Both the relative standard deviation 23 as the variation along the central 
pixel line^ 1 of T BGB are used to describe the global uniformity: 



SD(T^ B )* ± max(T*)-min(T* 



(11) 

central line 



(T* GB >* mean(T*) 

The variation of T as function of the A-position is estimated by the maximal standard 
deviation of T BGB of the blank films and the uniformly irradiated film. 

Cov (T RGB , A) = max (SD (T* GB )*) (12) 

Fllm S7,8,9,ll 

III. RESULTS 

A. Transmittance of blankfilms 

A blank film scan was selected for each month from January 13, 2010 until January 12, 
2011. For these films, major changes of both T* and the global uniformity are noticed, as 
shown in figure El The red and green channel (T BG )* ranges from 0.77 to 0.86, and from 
0.74 to 0.82. The blue channel has a larger range from 0.31 to 0.50. The global uniformity 
of these film ranges from 2% to almost 10%. The last two months an improvement of the 
blue channel's global uniformity is seen. 

B. Association between transmittance and delivered dose 



In a first exploration of the calibration curve, equation la, (3^ is set to zero, see figure [3j 



All the segments of Film^.^ are pooled for the /3jR GB -estimation. For the red channel all 
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parameters have highly significant results (p < 0.01); (5\ = 1.833 ±0.069, f3 2 = 0.199 ±0.013, 
and /?3 = 2.159 ± 0.083. The residuals show a slight dependence on the X-location on the 
film. Additionally a larger spread of the residuals for X = —9cm is noticed. 

The green channel has similar results with parameters /3i = 4.487 ± 0.340, 02 — 0.084 ± 
0.040, and 3 = 5.431 ±0.414 with significance levels p < 2-10~ 16 , p < 0.05, andp < 2-10~ 16 . 
The green channel residuals have a larger spread for the X = —9cm segments, but the X- 
dependence as shown in figure [3] is not present. 

The model was not appropriate for fitting the blue channel data, the /^-estimations are 
not significant. 

C. Storage conditions 

The environmental condition of our storage room is stable with an average relative humid- 
ity of 35.0±4.9%. The average temperature is 23.1±0.2°C According to the manufacturer 
the vacuum packing environment had the following stable conditions; relative humidity = 
35 — 55%, and temperature = 18 — 23°C. 

The evolution of (T*)* of the vacuum packed films is illustrated in figure [5] There is an 
increase from Film^ to a maximal value for Filniy 7 , which is the first film covered with an 
interleaving tissue. These maximum deviations are 1.25, 1.45, and 1.66% for the red, green, 
and blue color channel respectively. The difference is clearly dependent on the location in 
the stack, and whether an interleaving tissue was used or not. When the films are stored 
separately, the impact of the location in the transportation stack and the effect of the 
interleaving tissue reduces. The differences decrease over time to values smaller than 0.5% 
for the red and green channel and smaller than 1% for the blue color channel on the fifth day. 
In absolute values, (T^)* increases from 0.79 on October 20, 2010 over 0.80 on October 21, 
2010 to 0.81 on October 25, 2010. The same evolutions are noticed for (Tq B )*, with values 
of 0.77, 0.78, and 0.80 (green channel), and 0.45, 0.46, and 0.47 (blue channel). 

D. Six or eight segments 

The residual errors, £j, are calculated for the different data sets, j. Only data-sets 
which are pooled from a single calibration film, or from two calibration films are considered. 
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The difference between the eight and six segment approach is illustrated in figure [6j using 
\\ e j\\ = / Yl £ lj- Avoiding the segments with x = —9cm reduces ||e||. On average 

\j isdata set j 

these reductions are 0.43 ± 0.15, 0.28 ± 0.38, and 1.63 ± 1.05Gy for the red, green, and blue 
channel. An exception for this reduction is found for the film with the opposite MU-order, 
Filing. The 33 and 66 MU-segments of Filing are located on the X = —9cm position 
(table [I]). Eliminating these low dose segments results in an ||e|| increase for the green and 
blue color channel. 



E. Calibration 

Only the six segment approach is discussed in this section. 

1- PiRGB m estimations 

Table [IT] summarizes the significant /3;RGB-estimations. For the data sets pooled from 
a single film, significant red channel /^-estimations were found for Filing and Filing, 
the films with an arbitrary mix of low and high dose values on the different locations. 
Additionally, the blank segment of Filni£ 10 resulted in a significant To = ^-estimation for the 
green color channel. The green T^ = /^-estimation was borderline significant (p = 0.0585). 

For the calibration sets pooled from two films, the /^-estimations become highly significant 
for the red color channel (p < 0.01). Additionally, a mix of low and high dose values on dif- 
ferent locations results in significant estimation of all four /3jG-parameters (FilniE 3 -|- Filing, 
p<0.05). 

For this data set, Filing +Filni£ 4 , the blue channel T could not be determined, /3 1B and 
/3 3B are not significant, this in contrast to X^ = /3 2 b- 

2. Calibration curves 

In figure [7]the calibration curves resulting from the data sets pooled from Fi\m.E 3 -E 5 ,andE 10 
are compared with the calibration curve of Filing +Filni£ 4 . The results are evaluated with 
a 0.02Gy threshold, 1% of 2Gy, a commonly prescribed dose. Large differences (>0.02Gy) 
between the curves are noticed. Film£ 10 has the best performance. Filing (increasing MU, 
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with increasing segment number) under estimates the dose for values larger than 2Gy and 
over estimates the dose between 0.1 and 1.4Gy. Filing (opposite MU-order of film Film^) 
has an overall underestimation of the dose (A <0Gy), the low dose range (values <0.5Gy) 
has the worst performance. Filing (arbitrary MU-order) is comparable with Film^ but over 
estimations are seen where FilniE 3 has underestimations and visa versa. Filni£ 10 (=Filme 4 , 
but with a blank segment and higher dose for segment 6) has the better performance for the 
whole dose range. 

F. Error analysis 

1. Transmittance accuracy 

Figures [8] and [9] show respectively the local and global uniformity. The local uniformity 
is smaller than 1.08%, 0.97%, and 1.15%, for the red, green and blue channel. The global 
uniformity ranges from 2% to 7%. The blue channel has worse global uniformity than the red 
and green channel. The transmittance covariations with the X-position (Cov (T, X)) are 
5.7- 10~ 3 , 5.1 • 10~ 3 , and 5.5 • 10 -3 (red,green,and blue channel) for the six segment approach. 
Compared to the eight segment approach, the six segment approach has a 33% reduction 
of the red channel Cov (T, X). The green and blue channel covariances of the different 
approaches are comparable, ot was not found to be different for the two approaches. 

2. Residual dose error and dosimetric accuracy (e and <Td) 



Figure 10 illustrates the dose values based on the calibration curve of the Filing +Film B4 - 
data set (six segment approach). The dose from the film measurements, -DFiim rgb, and the 
corresponding IC measurements, -Die, are shown in the left figure. On the right side, the 
residual errors are displayed. An over estimation is found in the green channel for the un- 
irradiated segments. In the [0.04,2.5]Gy dose range, excluding all unirradiated segments, 
the red channel residual errors deviate 0.62±1.79% (meanilSD) from the expected dose, 
Djc, with a range of [—2.3%, 4.9%]. In the same dose range the green channel residual errors 
deviate 2.1±3.63% from the expected dose. These residual errors have a [—3.6%, 6.9%] range, 
with three outliers of 7.6, 9.53, and 16.6%. 

The dosimetric accuracy, an, is calculated for X = 0cm. For the red channel the calcu- 
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lated dosimetric accuracy is a good estimate for the upper limit of the residual errors. For 
the green channel the calculated ctd is an under estimation of this upper limit. 

IV. DISCUSSION 

A. T 

1. Stability 

Figure [4] illustrates an instability for the gafchromic EBT 2-films used in our department. 
The TJ -values are largely dependent on the production (variation between different batches). 
The instability is present in all three color channels but it has different effects, it stresses 
the necessity to calibrate each batch. An advantage of our calibration approach is the fact 
that the estimated calibration parameters indicate the zero dose transmittance (To), and 
therefore, inform the users about the changes of the film. 

For the same films, global non-uniformities were noticed. To work with a classic optical 
density protocol, a location dependent 0*-estimation should be made for each bach (f.e. 
Garcia-Garduno^ et al. proposes an average scan of five randomly selected blank films as 
0*-estimation for gafchromic EBT films). With a throughput of more than one box of films 
each month our protocol needs to be as transparent as possible. Therefore, we do not prefer 
the cumbersome 0*-estimation, and T is used instead of OD. 

2. Storage 

A stack of films without interleaving tissues allows the films to slide over each other, 
and static electric charges are created on the film. The resulting electrostatic forces stick 
the films to each other. This way the films are only exposed to the environment at the film 
edges. The interleaving tissue, avoids the direct contact of the films and serves as a guide for 
the humidity (David Lewis, ISP^, private communication). This is shown in figure pi where 
the films separated by interleaving tissues have similar To-values, and the films without 
interleaving tissue have increasing To-values in function of the location in the stack. When 
exposed to our local environmental conditions the films tend to reach a stable state, relative 
to each other, after 1 day. Meanwhile, the manufacturer reintroduced the interleaving tissue 
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for the transport of the films. 

B. Calibration curve 

1 . Protocol 

Only a correlation between e and the X-location was noticed (figure wl). The absence of a 
correlation between e and Die validates the assumption of using only eight IC measurements 
to describe the low dose contributions. 

The low local uniformity values (figure[8]) support the choice of the (small) 4x4 calibration 
fields. On the other hand, because of the absence of a o"r-difference between the six and the 
eight segment approach, ot is not affected by the global non-uniformities. Therefore, the 
lxl cm segment are small enough to evaluate local uniformity. 

To cope with global non uniformities a location dependent extension is added to the 
rational function, additionally the six segment approach was introduced. This approach 
decreases the degrees of freedom (six vs eight data points per film for four fa RGB-parameters). 
Nevertheless, this approach is an improvement for the protocol, indicated by the reduction 
of ||e||, figure [6j 

Figure [7] illustrates the better performance of the calibration curve from Filing. The 
improved significance of /3jG _ P arame ters of Filing compared to Filing indicate an improve- 
ment of the protocol due to the introduction of a no dose segment. The necessity of low dose 
values is further supported by the Film^-results, where the absence of low dose data points 
in the six segment approach enlarges the residual error, ||e||. This puts the dose levels used 
in the protocol in question. Future improvements of our protocol are the introduction of 
MU segment, and the other MU values will be selected to have equal intervals between the 
(T£)* data points. 

2. Location dependence 

Remarkable are the significance levels for the /3iRGB-parameters. For the single film 
calibrations (Film.E 1 -E (i ,E 10 ) , only the films with an arbitrary spread of the MUs on different 
locations generate significant results for the red channel. 
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For the two film methodology, the Filing + Film^-data set, results in significant /%- 
estimations for both the red and the green color channel. This in contrast to the Filni£; 2 + 
Film^-data set, where only the red /^-estimations are significant. The difference between 
these two data sets is the MU-order used to irradiate the films. Film^-I- Film^ has an 
opposite MU-order on the two films while Filing + Filing has the same MU-order on the 
two films. Therefore, the spread of high and low dose segments on different positions is 
found to be relevant for the calibration of the gafchromic EBT 2-films. 

In other words, the calibration curve is location dependent (D = f (T,X,Y)), which is 
also mentioned by Micke, Lewis, and Yii™. The study of such an effect, combining the 



dose and the X- location (e.g. a /3s A ■ D-term), was outside the scope of is this work. 
However, using two films with a spread of different dose levels on different locations, makes 
our calibration curve an acceptable surrogate for such a location dependent calibration curve. 

C. Error analysis 

1. Transmittance accuracy 

Although, T is used instead of OD, it is interesting to compare the uniformity results 
with literature values. After all, both T and OD form the input of a calibration curve, and 
the error on these input will determine the error of the resulting dose values. 

Excellent local uniformity values were found (< 1.15%) which is lower than 3% and 5% 
relative standard deviations for the MD-55-2-film listed in the TG 55 report?^. Zhu et 
a i[23] re p 0r t s 7%-15% global OD- variations for the MD-55 film, which is comparable with 
our transmittance results (±7%). Nevertheless, adaptations were necessary to cope with the 
global non uniformity (six segment approach, and A-dependent parameter in the calibration 
curve) . 

2. Residual dose error and dosimetric accuracy (e and ap) 

In a valid calibration protocol e should have a normal distribution with a zero mean 
and a standard deviation from the previous section (ed ~ N(0, <7d)). This relation is not 
evaluated because of the small amount of data, n = 6 ~ 8 per dose level. Therefore, we 
have chosen an over estimation of cr D , introduced by using the maximum of the local and 
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global uniformity variations in equations 10 and 12 

The red color channel cr^-prediction was able to estimate the maximal dose errors, figure 

For the green channel the calculated cr^- values are not sufficient to estimate the maximal 
dose errors in the protocol. Neither was the protocol capable to convert Tq q to dose, see 



figure 10 To cope with this problem, blank segments could be introduced in the calibration 
protocol. The significance of the (5\q and /^-estimations for the film with a blank segment 
(Filni£ 10 ) supports this assumption. 

V. CONCLUSION 

The proposed calibration protocol for absolute gafchromic film dosimetry has a straight- 



forward association between transmittance and dose, equation |2aj The approachable proto- 
col required only eight preparation measurements, and a single reference measurements on 
the day of calibration. The use of simple static 4x4 fields is suited for this purpose even as 
the use of IC-measurements to map the low dose contributions. 

Film changes are reported both, due to the production process as due to the storage 
environment. The strength of our protocol is that each calibration characterizes the physical 
parameters of the films, To, T^, and a factor scaling the impact of the dose ( J-). 

However, because of non-uniformities the original intended association between T and D 
required an adaptation. An additional term (/3 4 • X) and the avoidance of certain segments, 
are positively evaluated as improvements of the protocol. 

The protocol requires a spread of low and high dose segments on two calibration films 
and the segments with a higher spread of the residual dose errors should be excluded. 
Additionally the use of a blank segments is advisable. 

All color channels are dose dependent, but a trustworthy calibration protocol with error 
prediction was only found for the red channel in the ]0, 2.5] Gy dose range. In this dose range 
the red channel dose error range equals [—2.3%, 4.9%]. For the red color channel an upper 
limit for the dose error could be predicted using an error propagation analysis. The green 
color channel has the higher residual errors, and requires more attention to the estimation of 
Tq. The blue color channel is dose dependent but could not be calibrated with the proposed 
protocol. 
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Summary of our current calibration protocol 

Preparation measurements Eight preparation measurements are performed once to es- 
timate the low dose contributions. 

Storage The films are stored in separate folders of a file cabinet to ensure the same 
environmental exposure off all films. 

Blank Check After a few days of storage, for stabilization purposes, all films are scanned 
unirradiated. Thereafter, an automated script calculates (Tr GB )* ± SD {T^ GB )* , as well as 
the local uniformity off all the films. 

Dose levels A blank segment is introduced for a better estimation of To. In our current 
protocol the dose range is extended from to ~ 4.5Gy , using MU ranging from to 393. 
Where the number of MUs are chosen to be equally distributed in the red transmittance 
domain. 

2 Films Two calibration films are irradiated using opposite MU-orders. 

Calibration Check Before using clinical films, the calibration films are scanned and eval- 
uated. 

Segment Selection Based on the Blank Check and the Calibration Check a decision is 
made wether certain segments should be avoided or not (e.g. segments with X = —9cm). 
In our latest clinical evaluations such a segment avoidance was not necessary (batch Nr. 
A09031001B, A11011001, A12171002B). 

48h post irradiation We wait 48 hours before scanning the films. 

Scan all films in one session All the films that need to be evaluated are scanned suc- 
cessively, including the calibration films. 

Convert to dose A calibration curve is created for the scan session and all the films are 
converted to dose. 

Calibrate each box of films 

Acknowledgements Mr. K. Poels is acknowledged for extensive discussions. The 
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work. 
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Appendix A: Transmittance vs dose 



Equation la (Al) can be transformed to lb (A2) using the following routine juggling 



r= ft + ft ; g + ft.jr 



fo + D 
(T - fa • X) • (#, + D) = fa + fa ■ D 

(T-P4'X)'fo + (T-p i -X)-D = p 1 + fa-D 



{T-fa-fa-X)-D = fa-fa-{T-fa-X) 

fa - fa • T + fa fa • X 

T-fa-fa-X 
fa-T-fafa-X- fa 



D 
D 



fa-T+fa-X 



(Al) 



(A2) 



These equations can be reduced further using three physical parameters. The first two 
parameters characterizes the film using the zero and infinite dose transmissions (T and 
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Too). The third parameter, 3 scales the impact of the dose. 

01 + 02-D 
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D 



h-T- 3 4 • X - gi 
2 -T + 4 -X 
4 



ft 



*+&■£ 



1 + 
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+ /3 4 -X 



£> = & 



t - g4 • x - | 

& - T + 04 ■ X 
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03 



lim 



0i + 02 ■ D 
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i#2 (l'Hopital's rule) 
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T = ^ ^- + /3 4 • X 



D 



1 + 



ZJ 



£> = /?; 



3 • 



T - fa • X - Tp 

Too - T + /3 4 • X 
4 



^^+"+^■1 



1 + 



ft 
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(Tp + A .X)-T 
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TABLES 



TABLE I. Summary of the different irradiated and blank films used in this work. 



Experiment 


Film 


Batch Nr. 


Description 


Storage 


Vx-n 

V7-15 


A081610 
A081610 


Vacuum packed films without interleaving tissue. 
Vacuum packed films with interleaving tissue. 


Calibration 


E X 

E2,3 

E* 
E 5 

Eq 

E7 

Es,s,n 

Eio 


F03161001 
F03161001 
F03161001 

F03161001 

F03161001 

F03161001 
F03161001 
F03161001 


Calibration film with: MU = [33, 66, f 00, 125, 150, 175, 200, 225] and A,*Cai 100 ml 

Calibration film with: MU = [33, 66, 100, 125, 150, 175, 200, 225] and A,*Cal 100 ml 

Calibration film with the opposite MU order: 

MU = [225, 200, 175, 150, 125, 100, 66, 33] andA.iCai 100 mu = 1.0689Gy. 

Calibration film with a random MU order and A.iCal 100 mu = 1.0689Gj/; 

segment 1 to 8: MU = [150, 33, 225, 66, 200, 100, 125, 175]). 

The film is irradiated with ± equal number of MU's per row 

A,iCai 100 mu = l-0689Gy; 

segment 1 to 8: MU = [33, 225, 200, 66, 100, 175, 150, 125]). 

The film is uniformly irradiated with a 30x30 field, 2.38Gy at the central axis 

These films are left blank. 

Calibration film with the opposite MU order of E% and 

A,iCal 100 mu = 1.0689Gy 

this is the same order as E4, but segment five was left blank, 

and segment six got 100 + 125MU; 

segment 1 to 8: MU = [225, 200, 175, 150, 0, 225, 66, 33]) 


Stability 


Bi-13 


variate 


Blank films used in clinical practice through out 2010, see figure l 


\ and L 


1 
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TABLE II. Summary of the /3j/jGB-estimations. 



Film 




Pi 
meaniSD 


P[%] 


meaniSD 


P[%] 


ft 

meaniSD 


P[%] 


/3 4 
meaniSD 


P[%] 




Red 


1.64±0.04 


(0.07) 


0.24±0.01 


'0.08) 


1.9310.06 


(0.08) 


(-10±1)-10~ 4 


(0.40) 


E 5 


Green 


4.16±1.43 


(10.0) 


0.13±0.15 


'49.5) 


5.0211.77 


(10.6) 


(- 3±3)-10~ 4 


(42.5) 




Blue 


2.30±0.77 


(9.59) 


0.19±0.05 


'5.90) 


5.6411.91 


(9.86) 


( -9±1)-10~ 4 


(0.96) 




Red 


2.11±0.09 


(0.18) 


0.15±0.02 


1.05) 


2.471011 


(0.18) 


(-10±1)-10~ 4 


(0.99) 


E±o 


Green 


3.07±0.52 


(2.77) 


0.24±0.06 


'5.85) 


3.67±0.63 


(2.81) 


( -1±3)-10~ 4 


(68.0) 




Blue 


2.11±1.25 


(23.3) 


0.19±0.09 


15.6) 


5.15±3.05 


(23.4) 


( -6±2)-10~ 4 


(10.5) 




Red 


2.06±0.22 


(< 10- 3 ) 


0.15±0.05 


1.21) 


2.41±0.26 


(< 10- 3 ) 


(-15±4)-10~ 4 


(0.60) 


E2 + E% 


Green 


3.42±0.82 


(0.31) 


0.19±0.11 


12.3) 


4.05±0.97 


(0.31) 


( -7±5)-10~ 4 


(22.2) 




Blue 


2.03±1.64 


(25.1) 


0.21±0.12 


12.5) 


5.01±4.06 


(25.2) 


( -7±4)-10~ 4 


(16.1) 




Red 


1.90±0.13 


(< 10- 5 ) 


0.19±0.02 (< 


; 10- 3 ) 


2.22±0.17 


(< 10- 4 ) 


(-llzbl)-lO- 4 


(< 10- 3 ) 


E$ + £4 


Green 


2.90±0.54 


(0.07) 


0.26±0.06 


'0.28) 


3.43±0.67 


(0.09) 


( -6±2)-10~ 4 


(2.77) 




Blue 


1.43±0.67 


(6.50) 


0.25±0.04 


;o.o4) 


3.5111.68 


(7.08) 


( -8±1)-10~ 4 


(0.06) 
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FIGURE CAPTIONS 



FIG. 1. Schematic overview of the calibration image. The outer rectangle represents the 
film. All length indication are in [cm]. 

FIG. 2. The transmittance (T BGB ) allonge the central pixel line, of a uniformly irradiated 
film (table [II Film^). A reduction of the measured signal is noticed for X < —5cm. 

FIG. 3. A fit of the red color channel data of all the films: FilmB 1 _^ 11 . On the left the data 
is shown together with the fit Tr = ^^iscj+D IC • On the right the residuals are shown in 
function of their X location on the film. Note the X-dependence of the residuals as well as 
the larger spread of the residuals of the segments with X = —9cm 

FIG. 4. Blank film stability and uniformity from January 13, 2010 until January 12, 
2011. The upper part contains the average transmittance, (7r GB )*. The global uniformity 
( Agb)* ± maX m2n(T'" (r * ) I central line J is illustrated in the lower part. 

FIG. 5. Effect of environmental conditions of our storage procedure on the vacuum packed 
films. Average transmittance values, (T BGB )*, are presented relatively to the average trans- 
mittance of Film^ on October 20, 2010. 

FIG. 6. Eight versus six segment approach. The square root of the squared sum of the 
residual errors of the six segment approach (||e 6 1|) is subtracted from that of the eight 
segment approach (||e 8 ||). 

FIG. 7. Comparison of different calibration curves for the red channel transmittance. The 
curve of the calibration data set pooled from Filing +Filni£ 4 is considered as the gold 
standard, this curve is subtracted from the curves from Fimi£ 3 _£ 5 , and Film^ 10 . The 1% 
dose error, compared to 2Gy a common prescribed dose, is illustrated in red. 
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FIG. 8. Local uniformity: the segments' average transmittance and corresponding stan- 
dard deviation ((T RG b) ± SD (T RGB ) ) are shown for the following films: a calibration film 
(FilniE 4 ), a uniformly irradiated film (Filing), and a blank film (Fihng; g ). 

FIG. 9. Global uniformity ( ,l t RG ?J ± ^^ — rr^\ — (central line I of a uniformly irradiated 

V V RGB' mcan U ' I 

film (Filing) and three blank films (Filmg 8911 ). 

FIG. 10. Left, calculated dose values of Film^.^ using calibration curve based on E^-\-EA 
(six segment approach). Right, dose differences (e) of the left figure and dose standard 



deviation (ctd) calculated using equations 3a till|3d 
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FIGURES 




FIG. 1. Schematic overview of the calibration image. The outer rectangle represents the film. All 
length indication are in [cm]. 
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FIG. 2. The transmittance (7r GB ) allonge the central pixel line, of a uniformly irradiated film 
(table III Filing). A reduction of the measured signal is noticed for X < —5cm. 




1 .5 2.0 2.5 



FIG. 3. A fit of the red color channel data of all the films: Filing _£j n . On the left the data is 
shown together with the fit Tr = ' 2 Xa +D IC ■ On the right the residuals are shown in function 
of their X location on the film. Note the X-dependence of the residuals as well as the larger spread 
of the residuals of the segments with X = —9cm 
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FIG. 4. Blank film stability and uniformity from January 13, 2010 until January 12, 
2011. The upper part contains the average transmittance, (T^ GB )* . The global uniformity 
— ' '"''''' ± max ^ „„/™™ — - 1 central line ) is illustrated in the lower part. 



\Trgb)* 



mcan(T*) 





1 


1 


1 




























n 


20/ 10/20 10 







?M 


J/2010 


























" 


+ 


25/ 


" -SJ 1 u 


* 






I 


+ 


X 


* 


* 


* 


* 


* 










+ 
































* 








+ 






o 














. i 


O 





ft 


I 


$ 


8 

n 







□ 





□ 


$ 


• 


• 


- 


- 








R 


□ 


n 




□ 


D 


u 


n 


□ 


n 


- 








N 
























. 


B 


Z. 



























FIG. 5. Effect of environmental conditions of our storage procedure on the vacuum packed films. 
Average transmittance values, (7r GB )*, are presented relatively to the average transmittance of 
Film^ on October 20, 2010. 
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E2 E3 E4 E5 E6 E7 E3+E4 E2+E3 



FIG. 6. Eight versus six segment approach. The square root of the squared sum of the residual 
errors of the six segment approach (||£6||) is subtracted from that of the eight segment approach 

(Nil). 
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FIG. 7. Comparison of different calibration curves for the red channel transmittance. The curve of 
the calibration data set pooled from Filing +FilmE 4 is considered as the gold standard, this curve 
is subtracted from the curves from Film£ 3 _£ 5 , and Filmg; 10 . The 1% dose error, compared to 2Gy 
a common prescribed dose, is illustrated in red. 
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FIG. 8. Local uniformity: the segments' average transmittance and corresponding standard de- 
viation ((Trqb) =t S'D(Trgb) ) are shown for the following films: a calibration film (Filmg 4 ), a 
uniformly irradiated film (Film,E 7 ), and a blank film (Filing). 
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FIG. 9. Global uniformity I , };, RG ?J ± ri ^ — ' ^P — (central line ) of a uniformly irradiated film 

V \ 1 RGB/ mean(i ) 

(Filni£; 7 ) and three blank films (FimiE 8g u ). 



34 




FIG. 10. Left, calculated dose values of FilrnE 1 _£ 11 using calibration curve based on E% + E4 (six 
segment approach). Right, dose differences (e) of the left figure and dose standard deviation (o\d) 



calculated using equations 3a till 3d 
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